On isotropic turbulence in the dark fluid universe 
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As first part of this w ork, ex perimental information about the decay of isotropic turbulence in 
ordinary hydrodynamics, u 2 (t) oc t~ 6 ^ 5 , is used as input in FRW equations in order to investigate 
how an initial fraction / of turbulent kinetic energy in the cosmic fluid influences the cosmological 
development in the late, quintessence/phantom, universe. First order perturbative theory to the first 
order in / is employed. It turns out that both in the Hubble factor, and in the energy density, the 
influence from the turbulence fades away at late times. The divergences in these quantities near the 
Big Rip behave essentially as in a non-turbulent fluid. However, for the scale factor, the turbulence 
modification turns out to diverge logarithmically. As second part of our work, we consider the 
full FRW equation in which the turbulent part of the dark energy is accounted for by a separate 
term. It is demonstrated that turbulence occurrence may change the future universe evolution due 
to dissipation of dark energy. For instance, phantom-dominated universe becomes asymptotically a 
de Sitter one in the future, thus avoiding the Big Rip singularity. 

PACS numbers: 98.80.-k,04.50.+h,11.10.Wx 



I. INTRODUCTION 

Consider a spatially flat FRW universe with H = d/a the Hubble parameter. In standard notation the FRW 
equation 

-=-\n 2 {p + Zp) (1) 
a 

with k 2 = 8ttG, implies that for the scale factor a(t) to depict a curve concave upwards when drawn as a function of 
t, i.e. a > 0, it is sufficient that the equation of state (EoS) satisfies the condition 

P<-\p- ( 2 ) 
In order to obtain H > it follows, however, from the equation 

H = -~K 2 (p+p) (3) 

that p has to satisfy the stronger condition 

p<-p. (4) 

This is the phantom region, corresponding to a positive tensile stress in the dark energy fluid. It is known that 
phantom-dominated universe usually enters to finite-time future singularity (Big Rip) (see Refs. (U-IH)- The region 
— p < p < — p/3 is called the quintessence region, where the expansion of the universe is accelerated a > but not 
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super-accelerated H > 0. Note that effective quintessence universe may also end up in one of three possible types of 
future singularity Q. 

This paper analyzes the possibility whether there can be a turbulent micro structure superimposed on the dark fluid. 
Quite obviously it is the isotropic version of turbulence theory which then becomes most relevant, as one wishes not 
to disturb the macroscopic isotropy. Turbulence generally implies that there is a loss of kinetic energy into heat. 
Whereas this loss is often described macroscopically, in terms of a bulk viscosity C (cf., for instance, Refs. [4-6]), our 
program here is to replace £ with a microscopic shear viscosity. In other words, we wish to replace the length scale 
associated with the macroscopic £ by the Kolmogorov microlength scale, conventionally denoted by r/. 

The next section gives a brief overview of isotropic classical turbulence theory, in order to put our present approach 
in a proper context. Thereafter, we consider how the turbulence effect can be taken into account in the cosmological 
formalism. We focus on the following two possibilities: 

1) The turbulence effect can be included by the addition of a constant fraction, called /, to the laminar ordinary 
energy density p in the first FRW equation. Assuming / to be a small quantity, a perturbative solution to the first 
order in / can conveniently be found. 

2) Our second approach is to write the total energy density as a sum of four different parts: (i) a laminar dark energy 
part, (ii) a turbulent dark energy part, (iii) a radiation part, and (iv) an ordinary matter part. The fate of the universe 
is in principle predictable, on the basis of the weight given to each of the constituents of the total energy density. 

Options 1) and 2) are considered in Sections IIIII and ITVl respectively. 

We mention finally that our method of decomposing the fluid into turbulent and non-turbulent parts is in principle 
similar to the method recently used by Balakin and Bochkarev Q- These authors divided the energy as well as the 
pressure of the cosmic fluid into two components, one component referring to dark matter, the other referring to dark 
energy. 



II. EXTRACTS FROM KOLMOGOROV'S ISOTROPIC TURBULENCE THEORY @, @ 

Turbulence generally implies a loss of kinetic energy into heat. As mentioned, we replace the macroscopic bulk 
viscosity with a microscopic shear viscosity, corresponding to the Kolmogorov length 



V=[~) • (5) 

Here v is the kinematic microscopic shear viscosity and e is the dissipation per unit time and unit mass. 

Let I denote the external scale of the turbulence, with l/l the corresponding wave number. The large eddies move 
around with only a little dissipation of energy. According to the so-called second hypothesis of Kolmogorov Q , in an 
isotropic region the motion is entirely determined by friction and inertia. There occurs a continuous flux of energy 
transferred by means of a hierarchy of eddies corresponding to the dissipation e. Let A characterize the size of an 
eddy, k = 1/A being the corresponding wave number. The equilibrium range is that for which all memory of the flow 
is lost, 

fc»l/Z. (6) 

If u\ is the typical velocity of an eddy of size A, the internal Reynolds number is Re a ~ Xu\/v. For increasing 
values of k, ReA decreases. Dissipation becomes important when Re> ~ 1. This is just the condition leading to the 
Kolmogorov length ([5]). 

If the Reynolds number of the flow as a whole is high, the wave numbers 1 [I and I/77 are widely separated, and 
there exists an inertial subrange characterized by 

y<fc<-, (7) 
I i] 

in which the fluid behaves like a non-viscous fluid. A famous formula for the spectral energy density, conventionally 
called E(k), in the inertial subrange is 

E{k) = ae 2 ' 3 k^' 3 , (8) 

where a « 1.5 is the Kolmogorov constant. 

For practical purposes the von Karman interpolation formula for E(k), linking the region of small k to the region of 
high k, is useful in order to calculate the total energy density E by integrating over all wave numbers (cf., for instance, 
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Ref. [1(|). We abstain from going into further detail here. Of interest for us here is, however, is the empirical decay 
law for isotropic turbulence (of course, under the assumption that it is left to itself; there are no external sources). 
Based on grid experiments in wind and water tunnels, it turns out that the mean kinetic energy ^u 2 (i) decays as 

u2(t) oct~ 6/5 ; (9) 

cf. Refs. [U, [IH, as well as the theoretical treatment in [l(|. We shall make use of this relationship in the following. 
Consider now the classical equation of motion for a viscous fluid: 

d t {p m Ui) + d k U ik = 0, (10) 

where p m is the mass density and Hit the momentum flux density tensor 

n,: fc = pS ik + pmUiUk - p{dkUi + diU k ), (11) 

p being the shear viscosity. Taking the mean of this equation, observing that u,; = in homogeneous and isotropic 
turbulence, we get 

2 

U xx = Uyy = U zz = p ef£ =p + -pturb) (I 2 ) 

where p c g is the effective which takes into account that the thermodynamical pressure is augmented by a term 
(2/3)ptmb associated with the turbulent energy density, 

Pturb = ^P,nWu~i = ^PmU 2 . (13) 

Thus Pturb designates a mean quantity. 



III. FRW EQUATIONS WITH A FRACTION / OF THE ENERGY AS TURBULENT ENERGY 

We will now consider this fluid in a cosmological setting. As usual in turbulence theory we may start by decomposing 
the fluid velocity Ui into a mean component Ui and a fluctuating component u[ , u, = Ui + u^. However, in a comoving 
reference frame U — 0, so that we can simply replace is! i with Uj. In order to keep the turbulent part of the cosmic 
fluid separate from the non-viscous (non-turbulent) part, we shall from now on endow non-viscous quantities with a 
subscript zero. Thus Eq. (TT2"j) is rewritten as 

2 

Peg = PO + ^Pturb, (14) 

with pturb = |pmi 2 as before. Note that the turbulent energy in the comoving frame is regarded as a nonrelativistic 
quantity. The total energy density can be written as 

P = PO + Pturb- (15) 

As mentioned above, we shall first look for a solution of the cosmological equations when a constant fraction / of the 
energy exists in the form of turbulent energy. As we are primarily interested in the dark energy epoch of the universe, 
we assume henceforth that the thermodynamical parameter w satisfies the inequality w < —1 (see Eq. (|21j) below). 
The initial instant for our considerations will be denoted by ti n . 

By assumption we can thus write the energy density at time t = t- m as 

p(t in ) = P0(tin)(l + /). (16) 

We see that / can be interpreted as the ratio between the turbulent energy and the total energy at t- ln , 

Pturb 



/ = 



Pa 



(17) 



As po includes the rest mass, this shows showing that the assumption / -C 1 is a plausible one. 
For t > t{ n we now require pturb to decay with time as 



Pturb « t~ 6 '\ (18) 
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in accordance with Eq. (J9j) for ordinary turbulence. We can write p(t) for t > t- m in the form 

p(t)=p (t)[l + fp 1 (t)], 

with 



6/5 



(19) 



(20) 



When t = tin, this agrees with Eq. (pi 

The equation of state (EoS) for the cosmic fluid is now to be introduced. This can be done in various ways. We 
choose write it in such a way that only non-turbulent quantities are involved, 



Po{t) = wpo(t), w : constant. 
The effective pressure will analogously to Eq. (fl"9)) be expanded as 

PeS (t) =W P0 (t)[l + f Pl {t)]. 

At t = t- m it follows that pi(iin) is actually a known quantity, 

, \ 2 
wpi(t in ) = -. 

In the same way we can expand the scale factor as 

o(t) = a (t)[l + f ai (t)} , 

and analogously for the Hubble factor 

H(t) = H (t)[l + fH 1 (t)]. 

The correction terms {p\, a%, H\} are all of zeroth order in /. 
From H = a/a we get at once 

di = H Hi, 

whereas the first FRW equation H 2 = \n 2 p yields the first order relationship 



H X = - Pl . 



The Hubble parameter H{t) thus satisfies the equation 

H(t) = H (t) 



6/5' 



We still need to determine H n (t). It can be found from the non-turbulent FRW equations 



do , 1 



a 2 u 2 

which, together with the conservation equation for energy, T 0,y ; „ = 0, yield 

Pa + 3H (p +p ) = 0. 

From Eq. ((29]), H = {V3/6)Kp /^/p Q , and as H = -H 2 + a/a we get 



(21) 
(22) 

(23) 

(24) 
(25) 

(26) 
(27) 

(28) 

(29) 
(30) 

(31) 

(32) 
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Here we have for convenience introduced the symbol 7, defined as 

7 = 1 + w. 

The solution of this equation is 

Ho{tin) 



H (t) 



Thus 



Correspondingly, we obtain 



H(t) = 



P(t) 



Ho(t- SD ) 



l + f7flo(tin)(*-*in) 



6/5 



n 2 [l + f7-Ho(*in)(*-*in)] 



1 + / 



t 



6/5' 



(33) 



(34) 



(35) 



(36) 



We can now draw the following important conclusion: both for the Hubble factor, and the energy density, the influence 
from the turbulence fades away when t 3> tin- To a good approximation the future singularity occurs at the same 
instant t — t s as if turbulence were absent, i.e. 



t s t[ n -\- 



3|7l#o(*m) 



Both H and p diverge at t = t s . Near t s , as t in /t s <C 1, 

Hq [ti, 



H(t) 



1 - t/ts 



t -> t. 



(37) 



(38) 



Consider next the correction a\ to the scale factor. From Eqs. and (|27p. 



61 = \hoPi , (40) 

from which we obtain by integration, setting x = t/t- ln — 1, 

1 ft/txn-i dx 1 

«'W = 2 ff "<'*/ + - <41) 

We need not calculate this integral in full, but note that it diverges logarithmically at x = t s /t in . Omitting multi- 
plicative factors, we write for the dominant part 

ai(t)~ln(l-t/t,), t->t,. (42) 

Thus the modification coming from turbulence is in this case itself turbulent. As the solution in the case / = is now 

a ° {t) [l + hH (t m )(t-t in )]V^ ' (43) 

it follows however from the expansion (|24j) that the divergence in ao(t) is much stronger than the turbulence modifi- 
cation. The dominant term near t = t s is thus 
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just as in the case / = 0. 

So far, we have made use of the first FRW equation only; we have not considered the pressure in the cosmic fluid. 
To deal with the pressure, we have to take into account the second FRW equation also, for instance in the form 



dt 



(pa 3 ) = -3H PeS a 3 



By expanding in the parameter / in the same way as above, we obtain to first order 

p\ + 3di 



Pi 



Pi-Hi- 



Inserting for p±, Hi, a±, and Hq we get 



Mt) = - 2 



14 1 
w 5 H()Wt 



Pi(*) 



(45) 



(46) 



(47) 



It is here to be observed that if one extrapolates this expression back in time, until the initial instant t m , the expression 
does not in general agree with the previous equation (f23|) . The reason for this is that our condition ([23} on the initial 
pressure makes the system mathematically over-determined. Of most physical interest is, however, the cosmic pressure 
at late times, Hq and t large, in which case the last term in Eq. (|47| fades away and we get 



6/5 



(48) 



When w lies between -1/3 and -1, i.e. in the quintessence region, the value of pi is actually higher than when w < — 1. 



IV. TURBULENT DARK ENERGY DENSITY COMPONENT AS A SEPARATE TERM IN THE FRW 

EQUATION 

We now leave the perturbative approach, and consider instead the first FRW equation together with the energy 
conservation equation when the total energy density is written as a sum of four different parts: first, a dark energy 
contribution consisting of a laminar part Pdark and a turbulent energy part pturb so that 

Pdark energy Pdark + Pturb! (49) 

secondly, a radiation part p ra d and an ordinary matter part Pmatter- The FRW equation thus reads 

3 

~nH = Pdark + Pturb + Prad + Pmatter ■ (50) 
K 

As in the previous section, we follow the development of the universe from the instant t m onwards, and we adopt the 
same empirical law for the time development of the turbulent energy density, 

/ * v 6/5 

Pturb = Pturb(tin) I — J ■ (51) 

The time derivative will be written in the form 



where 



Pturb = -CplHh > ( 52 ) 



C= [pturbtUr 576 • (53) 



Now consider the energy balance equation in which ptmb is considered as a source term, 



Pturb + 3if (pturb + Pturb) = -CpJurb ' ( 54 ) 
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Here, we assume Eq. (|52|) should hold only in the flat universe. In the FRW universe, Eq. (|52|) should be changed as 
in Because of the turbulence, the kinetic energy changes into heat and the heat then becomes radiation. Then 

we may consider the the conservation law for radiation in the form 

Prad + 3H(p lad + p rad ) = Cp\Hl . (55) 

Here, p ra d = p ra d/3. For definiteness let us consider the case where the turbulent part pturb dominates. Then Eq. (|50[) 
gives 

We shall now assume that the EoS parameter Wturb of the turbulent part is constant. The EoS for the turbulent 
quantities is written in conventional form, 

Pturb = "^turb Pturb ■ (57) 

Then, with the definition 7t U rb = 1 + ttfturb we can write Eq. (|54l) as 

= Pturb + ^V3 7turb Ptufb + C ^turb ' ( 58 ) 

We shall discuss three different options for this equation: 

(i) If 

Pturb>(^) 3 , (59) 

the third term dominates compared with the second, and we recover the expression (|52[) . It means that pturb behaves 
as in flat spacetime, pturb oc i~ 6 / 5 . 

(ii) By contrast, if 



Pturb < 



m 3 ■ « 



the second term dominates compared with the third, and the turbulent term becomes negligible. Then pturb behaves 

as a usual perfect fluid, giving pturb oc et -37turb and H ~ 3 ^ 2 - j (it is here assumed that u> tur b > — !)• 

(iii)If 

Wturb < "I , (61) 

there exists remarkably enough a solution where pturb is a constant, 

/ (j \- 6 / 5 

Pturb= ■ (62) 

V OTturb / 

We now consider the case where the scale factor a is given as a function of the cosmological time, a = a{t). Equation 
(|54p may be rewritten as 

j t (a 3 ^ Ptulh ) = -CaT**-*' 2 (a 3 ^»Vtu lb ) ^ , (63) 

which can be integrated to yield 

-6/5 



Pturb(t) = a(t)- 3 ^ jf dt'a{t')-^' 2 + C ) 



(64) 



Here Co is a constant of integration, which can be determined from the initial condition at t = t ln . When a(t) is a 
constant, Eq. (|64| reproduces the standard result: pturb oc when t is large enough. In case of lUturb = — 1, even if 
a(t) is not a constant, we obtain pturb oc t~~s for large t, again. 
We can also integrate (|55|) to obtain 



t 



PradW = Ca{t)-' 1 ( j dt'a(t>y Pturh (t) +Cx). (Or.) 
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Here C± is a constant of integration. 

When Wturb 7^ —1, if we consider the case of the de Sitter space: a(t) = aoe Hot with constants and H , we obtain 

p turb = a"**"*** ( C U-***Sl*» _ e -**^t\ + Koa?±Z | . (66) 

I 37turb V I 07turb 

When w > — 1 and t is large enough, we find 

Pturb « e- 3 ^» bff «* . (67) 
On the other hand, when w < ~ 1 and t ^> tin, Pturb goes to a constant 

_ 6 

Pturb -► f-^~ ) ' , (68) 
V 07tmb / 

which corresponds to (|62|). 

When Wturb 7^ — 1, if we consider the case of an effective quintessence- like power law expansion, a(t) = aot ha with 
constants ao and ho, we find 

If we consider the case of the phantom-like power law expansion, a(t) = a (t s — t) °, we find 

Pturb — [t a — t) 

_ 6 

* ( o (1 + }L/,) - C - «' +IW ") ■ 0») 

„ _ _G _ 

Especially when ho = — 3 ^ - with lOturb < — 1) we find ptmb oc (t s — t) 5 when i — > t s . In this case, from (|65l) . 
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we also obtain p rac j oc (t s — t) 5 . Then both pturb and p rac j increase rather rapidly, although the rate of increase is 
smaller than for the energy density of the phantom dark energy where p p hantom oc (f s — t) . 

Just for simplicity, we now consider the turbulence of the dark energy with u>t U rb = — 1. Then Eq. f|64[) gives 



Ptnrb = (Jtt ) (' -In) . (71) 

Here to = t- m — and we assume t- m > to . We now consider the simple case where the contribution from the matter 
(except radiation generated by the turbulence) can be neglected and the non-turbulent part of the dark energy has 
the constant EoS parameter u>dark = — 1, as for the cosmological constant. We write pdark = A. 
By using (|7"Tj) . we may rewrite (|55|) as 

! fa 4 *-*) = (!) T <r»(t-td)-*« 4 , (72) 

since p ra d = p ra d/3. Then by multiplying the FRW equation ([50)) (pdark = A and p ma tter = 0) with a 4 and differenti- 
ating with respect to the cosmological time t, we obtain 

, 6 

6 \ 5 



-(i7 + 2^)^4A + 4(-j J (t-tort. (73) 

The second term decreases with time and therefore for large t, we obtain the asymptotic de Sitter universe, where H 
is a constant 



// =U Ll ^J±f. (74) 
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Let assume the turbulence begins at t = to. Then since we assume Wturb = w<jark = — 1, the total dark energy behaves 
as a cosmological constant and the de Sitter universe is realized, where 



H = H 



TO 



\ 



y(A+(^l -'-) r I IT-.] 



After the turbulence begins, the cosmological constant decays and H becomes smaller and at the late time, the 
universe reaches the asymptotic de Sitter universe with H = Hlq. Then we may assume Hiq 3> Hlq. 

Thus, the situation considered is one where there is matter with vanishing EoS parameter and phantom dark energy. 
Without turbulence there could occur a phantom crossing. After such a crossing, the dark energy dominates. If there 
is no non-turbulent part, the density of dark energy is large and finally satisfies the condition (|59[) . the dark energy 
dissipates and converts into radiation. Hence, the accelerated expansion will terminate. Or before satisfying the 
condition (|59p. the energy density goes to a constant (|62p. corresponding to the asymptotic de Sitter space. 

Some remark is in order. Let us imagine the inflation ended by the turbulence and the turbulent part pturb of the 
energy-density generated the reheating. p t urb could been converted to the radiation. Then the energy- density p ra d of 
the radiation after the reheating could be of almost the same order with pturb- Therefore, we may obtain 

6 

/ 6 V 

Prad ~ Pturb = I J (*in ~ *o) • (76) 

Here, Eq. (ITTj) is used. Applying the Stefan-Boltzmann law 

Prad = VT 4 , (77) 

with the constant a, one may evaluate the reheating temperature T as 

which may give a constraint on the parameters. This may indicate on observational manifestations of turbulence. 



V. DISCUSSION 



In summary, we have investigated the role which may be played by isotropic turbulence in a dark fluid universe 
at late times. Using experimental information about decay of isotropic turbulence in classical hydrodynamics, two 
different approaches to modification of FRW equations are proposed. The consequences of turbulence presence at 
future phantom/quintessence universe are studied. It is demonstrated that it may change the characteristics details of 
finite-time future singularity, for instance, the behavior of scale factor at Big Rip. However, it seems that turbulence 
cannot remove the future finite-time singularity in a perturbative approach. In a non-perturbative approach, the 
turbulence presence may change the evolution of dark fluid via its dissipation. In particular, it may terminate the 
accelerated evolution or convert the phantom-dominated universe into a future de Sitter space. Hence, inclusion of 
turbulence may suggest a way to resolve the future singularity problem. 

In principle, one should not limit oneself only to the known decay law for isotropic turbulence. It is quite possible 
that at large scales more general laws should be implemented. In this case, the description is somewhat similar to the 
inhomogeneous equation of state of the universe introduced in Refs. |LML5j |. Moreover, the explicit scenario of the 
(effective) turbulence emergence should be developed. Perhaps, the easiest way to realize it is to use of cosmological 
reconstruction in modified gravity (for a recent review, see [16J). In this case, the use of turbulence turns out to be an 
effective description due to a corresponding modification of gravity. Alternatively, the origin of turbulence in a dark 
fluid universe may be related to an inhomogeneous equation of state for the universe. 

The important question in connection with our proposal of including turbulence in a dark energy universe is of 
course related to the search of observational evidence for turbulence signatures. This is expected to be possible in the 
near future when observational data will give us more accurate information about the cosmological equation of state. 
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